UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADP013638 

TITLE:  Toward  the  De-Mystification  of  LES 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  DNS/LES  Progress  and  Challenges.  Proceedings  of  the  Third 
AFOSR  International  Conference  on  DNS/LES 

To  order  the  complete  compilation  report,  use:  ADA412801 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADP013620  thru  ADPO 13707 


UNCLASSIFIED 


TOWARD  THE  DE-MYSTIFICATION  OF  LES 


C.D.  PRUETT 

Department  of  Mathematics  & Statistics 
MSC  7803 

James  Madison  University 
Harrisonburg,  VA  22807 


1.  Introduction 

Some  35  years  have  elapsed  since  large-eddy  simulation  (LES)  was  intro- 
duced as  a computational  tool  for  weather  modeling  (Lilly,  1966).  Because 
direct  numerical  simulation  (DNS)  will  remain  prohibitively  expensive  well 
into  the  foreseeable  future  for  high-Reynolds-number  flow,  LES  continues 
to  hold  great  promise  for  simulating  flows  of  engineering  interest.  However, 
despite  the  attention  the  method  has  received  over  the  past  decades,  and 
despite  important  recent  developments  (e.g.,  the  notions  of  “explicit”  fil- 
tering and  localized  or  “dynamic”  modeling  (Germano  et  al.,  1992)),  LES 
has  been  slow  to  mature  as  a predictive  tool.  In  contrast,  parabolized  sta- 
bility equation  (PSE)  methodology,  for  example,  which  originated  in  the 
mid  1980’s  (Bertolotti  et  al.,  1992),  matured  quickly  and  is  now  ready  for 
use  by  the  aerospace  industry  for  transition  prediction. 

In  the  author’s  view,  the  relatively  slow  adaptation  of  LES  as  a pre- 
dictive tool  arises  not  from  any  fundamental  flaw  in  the  idea  itself  but 
most  likely  from  misconceptions  that  widely  permeate  the  practice  of  LES. 
The  current  paper  addresses  three  pervasive  misconceptions,  each  of  which 
reveals  a lack  of  clarity  regarding  the  properties  of  digital  filters  and  the 
relationship  between  the  grid  filter  and  the  subgrid-scale  (SGS)  model. 

2.  Misconception  1:  Grid  Filtering  is  a Mere  Formality 

From  an  otherwise  excellent  recent  text  on  turbulent  flows  (Mathieu  and 
Scott,  2000)  comes  the  following  quote: 

Note  that,  if  one  forgets  that  the  Ui  is  supposed  to  represent  the  filtered  velocity 
field,  the  filter  plays  absolutely  no  role  (emphasis  added)  in  the  final  LES 
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equations,  which  can  be  understood  in  terms  of  adding  a subgrid  viscosity  term 
to  the  Navier-Stokes  equations,  without  introducing  the  notion  of  filtering. 

The  mistake  here  is  confuse  formal  invariance  under  filtering  operations 
with  quantitative  invariance.  To  be  precise,  the  grid-filtered  Navier-Stokes 
equations  are  given  as  follows: 

duk  d(ukui)  _ 
dt  dxi 

Tkl 


op  , 1 d2uk  drki 


+ 


+ 


dxk  Re  dxidxi  dxi 

UkUl  - UkUl 


(1) 

(2) 


Whereas  the  form  of  Eq.  1 is  independent  of  the  filter,  the  value  of  the 
residual  stress  rki  is  not. 

To  illustrate,  suppose  that  grid  filtering  is  applied  at  a wavenumber 
cut-off  kc  greater  than  the  Kolmogorov  wavenumber  kn  (Fig.  1).  Then,  in 
effect,  one  conducts  fully  resolved  DNS  rather  than  LES,  in  which  case,  for 
all  practical  purposes  Tki  = 0.  As  written,  Eq.  2 is  missing  a parameter.  In 
reality,  rki  = Tki(kc).  That  is,  the  exact  residual  stress  depends  first  and 
foremost  upon  the  cut-off  of  the  grid  filter.  What  else  does  it  depend  upon? 
Figure  2 compares  exact  residual  stresses  computed  from  a priori  analy- 
ses of  isotropic  turbulence  (Pruett  and  Adams,  2000)  by  filtering  with  two 
different  grid  filters:  one  of  2nd-order,  the  other  of  4th-order.  The  correla- 
tion between  the  exact  residual  stresses,  which  differ  in  both  distribution 
and  amplitude,  is  only  0.5!  Thus,  at  minimum,  rki  depends  upon  both  the 
cut-off  and  the  order  of  the  grid  filter.  Can  the  dependence  of  the  residual 
stress  upon  the  filter  be  fully  characterized? 

Yes.  A filter  is  completely  characterized  in  terms  of  its  transfer  func- 
tion H{0,  which  expresses  the  action  of  the  filter  on  the  Fourier  harmonic 
exp (ikx),  where  £ = kAx.  That  is,  if  the  filter  is  represented  by  the  param- 
eterized (A)  convolution  integral  u{x)  = G(x'  — x,  A)u(x')dx',  then  H 
is  the  Fourier  transform  of  the  kernel  G,  and  G is  its  inverse  transform.  For 
example,  consider  the  top-hat  filter  whose  kernel  is 


G(x,  A)  = j 


1/A  if  |ar|  < A/2 
0 otherwise 


(3) 


Trapezoidal-rule  quadrature  over  a filter  width  of  A = 2 Ax  results  in  the 
discrete  3-point  top-hat  filter  of  weights  [1/4, 1/2, 1/4],  whose  transfer  func- 
tion is 


H(0  = 


1 + cos(£)  x2  x 4 

2 “ ~2(2!)+2(4!) 


(4) 


Alternately,  the  action  of  a filter  can  be  examined  by  Taylor-series  expan- 
sion of  a filtered  field  in  terms  of  the  un-filtered  field  and  its  derivatives. 
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For  example, 

u(x)  = u(x ) + a2ii"(x)Ax2  + <24  (a?)  Ax4  + ...  (5) 


(Because  we  have  presumed  G to  be  symmetric,  the  expansion  has  only 
even-ordered  terms.)  For  convenience,  the  factorials  in  the  expansion  have 
been  absorbed  into  the  coefficients  ft2m.  A filter  is  said  to  be  of  order  2m  if 
0.2m  (no  > 0)  is  the  first  non- vanishing  coefficient.  Equations  4 and  5 carry 
the  same  information  in  different  guises,  which  becomes  apparent  if  Eq.  5 
is  applied  to  the  function  u(x)  = exp (ikx),  from  whence  it  follows 


02m 


(2m)! 


(6) 


Thus,  the  coefficients  of  Eq.  5 are  completely  determined  by  derivatives  of 
the  filter’s  transfer  function.  For  the  three-point  top-hat  filter,  for  example, 

°2m  = 2(2  m)\- 

Now  let  Gx  and  Gy  be  discrete  filters  in  the  x and  y dimensions,  whose 
Taylor  coefficients  are  a2?n  and  62m,  respectively.  For  two-dimensional  fil- 
tering, suppose  u = Gy(Gx((u)).  By  two-fold  expansion,  one  obtains  the 
following  expression  for  the  exact  residual  stress: 

duk  dui  2 duk  dui  2 
-tm  = 2a2——Ax*  + 262 -^7- ^77  Ay 


+ 
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dy2  dy2 
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d3uk  dui 
dy 3 dy 


Ay4 


+0(he)  where  h € {A.t,  Ay} 


Eq.  7 is  completely  general;  it  is  valid  for  any  grid  filters  Gx  and  Gy  whose 
coefficients  are  known,  provided  the  associated  Taylor  series  converge.  The 
issue  of  convergence  is  addressed  in  Vasilyev  et  al.  (1998)  and  Pruett  et  al. 
(2001).  In  short,  convergence  is  guaranteed  for  positive  symmetric  filters. 

We  conclude  that  the  exact  SGS-stress  tensor  is  completely  determined 
by  the  grid  filter,  as  expressed  through  derivatives  of  its  transfer  function. 


3.  Misconception  2:  Any  Model  Goes  with  Any  Filter 

In  the  conventional  practice  of  LES,  the  grid  filter  and  the  SGS-model 
are  often  selected  independently,  as  observed  by  Piomelli  et  al.  (1988).  In 
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contrast,  previous  discussion  reveals  the  filter,  the  exact  SGS  stresses,  and 
the  model  to  be  fundamentally  interrelated. 

To  be  specific,  Fig.  3 compares  the  transfer  functions  of  filters  of  selected 
orders.  The  order  of  a filter  is  related  to  the  flatness  (i.e.,  the  number  of 
vanishing  derivatives)  of  its  transfer  function  at  £ = 0.  Spectral  filters, 
characterized  by  a sharp  cut-off  in  Fourier  space,  act  with  infinite  order. 
From  Eq.  7,  it  follows  that  the  leading  term  of  tm  is  of  order  h2rn  if  the  grid 
filter  itself  is  of  order  2m  (m  > 1). 

The  most  common  models  for  LES  are  of  eddy- viscosity  type  (Mathieu 
and  Scott,  2000),  of  which  the  Smagorinsky  model  (Ml)  is  the  most  famil- 
iar. Because  Ml  is  scaled  by  h2,  it  is  consistent  only  with  filters  of  (first 
or)  second  order.  In  particular,  Ml  is  inconsistent  with  spectral  filters. 
Even  when  used  in  conjunction  with  second-order  filtering,  the  Smagorin- 
sky model  correlates  poorly  with  exact  (E)  residual  stresses,  as  observed 
in  experiments  by  Liu  et  al.  (1994)  and  computational  studies  by  Pruett 
and  Adams  (2000).  Typically  C(E,M  1)  < 0.2,  as  suggested  in  Fig.  4.  In 
contrast,  similarity  (M2)  and  gradient  (M3)  models  perform  considerably 
better  in  a priori  analyses.  Typically,  C(E , M2)  > 0.8  and  C(E,  M3)  « 0.6, 
respectively  (Liu  et  al.,  1994;  Pruett  and  Adams,  2000). 

The  many  well-known  deficiencies  of  the  Smagorinsky  model  help  to 
explain,  in  physical  terms,  its  poor  performance.  Among  these,  Ml  is 
overly  dissipative,  which  undermines  its  successful  application  to  laminar- 
turbulent  transition,  and  it  is  isotropic,  which  limits  its  applicability  in  the 
near-wall  region.  However,  from  a mathematical  perspective,  the  deficiency 
of  Ml  is  apparent:  formally  it  doesn’t  match  the  leading-order  terms  of 
Eq.  7.  In  contrast,  M3  matches  at  leading  order  term  and  M2  matches 
at  leading  order  and  partially  matches  beyond  that.  The  author  is  led  to 
concur  with  Leonard  (1997)  that  the  Smagorinsky  model  has  “little  justi- 
fication” beyond  some  nice  properties. 

To  conclude,  further  advancement  of  LES  requires  that  more  attention 
be  devoted  to  filter-model  consistency,  to  models  that  have  greater  fidelity 
to  Eq.  7 (for  example,  similarity  models)  than  does  the  Smagorinsky  model, 
and/or  to  promising  new  mathematical  techniques  such  as  deconvolution 
(Domaradzki  and  Saiki,  1997;  Stolz  and  Adams,  1999)  that  avoid  explicit 
modeling. 

4.  Misconception  3:  Fixed- Width  Filters  are  Fine  for  LES 

Whereas  the  wavenumber  k*  and  the  grid  increment  Ax*  are  dimensional 
quantities  with  dimensions  [1/L]  and  [L],  respectively,  their  product  £ = 
k*Ax*  is  dimensionless.  (Here,  asterisks  denote  dimensional  quantities.) 
Several  problems  in  LES  result  from  a failure  to  non-dimensionalize. 
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To  be  specific,  the  maximum  wavenumber  k*mix  in  a simulated  flow 
is  determined  by  physical  considerations.  For  fully  resolved  DNS,  A:^ax  > 
k*  where  k*,  the  Kolmogorov  wavenumber,  depends  solely  on  Reynolds 
number.  For  LES,  k*nax  « k* . where  k*  lies  in  the  inertial  range  of  the 
energy  spectrum  (Fig.  1).  Either  way,  physics  determines  k^ ax. 

On  the  other  hand,  the  grid  increment  Ax*  is  mandated  by  numerical 
considerations,  that  is,  by  the  resolution  necessary  to  resolve  the  smallest 
eddies  for  the  numerical  scheme  of  choice.  For  example,  for  Fourier  and 
Chebyshev  spectral,  and  fourth-  and  second-order  finite-difference  (FD) 
schemes,  respectively, 

Ax*  = tt / k*nax  (Fourier  spectral) 

Ax*  = 2 / A.'*nax  (Chebyshev  spectral)  (8) 

Ax*  = 7r/(6fc^ax)  (4th-order  FD) 

Ax*  = 7r/(l 6&£,ax)  (2nd-order  FD) 

For  the  spectral  schemes,  the  stated  resolution  results  from  application  of 
the  Nyquist  criterion.  For  the  FD  schemes,  resolution  is  based  upon  the 
author’s  computational  experience  (Pruett  et  al. , 1995). 

A fixed-width  filter  (e.g.,  the  3-point  top-hat  filter)  has  the  disadvan- 
tage that  its  cut-off  is  “hardwired”  to  the  grid  increment  Ax*.  In  contrast, 
a tunable  filter  permits  the  cut-off  k*  to  be  adjusted  independently  of  Ax*. 
Figure  6 presents  the  transfer  functions  of  a family  of  one-parameter  second- 
order  filters  of  Pade  type  (Lele,  1992),  as  functions  of  their  dimensionless 
cut-off  £c  = k* Ax*,  where,  by  definition,  H(£c)  = 1/2.  The  filter  is  con- 
tinuously adjustable  over  its  entire  domain  of  cut-off  values  0 < £c  < n, 
which  represent  varying  levels  of  numerical  dissipation.  In  particular,  fc  = w 
provides  no  dissipation,  and  « 0 provides  maximum  dissipation.  The 
fixed- width  top-hat  filter  corresponds  to  the  particular  value  £c  = it/ 2. 

What  then  are  the  problems  with  fixed-width  filters?  First,  for  FD  nu- 
merical schemes,  allowing  the  dissipation  of  the  filter  to  be  set  by  the  grid 
increment  virtually  guarantees  that  the  unintended  truncation  error  of  the 
numerical  scheme  will  contaminate  the  intended  dissipation  of  the  filter. 

Second,  suppose  we  wish  to  compare  LES  against  the  results  of  a priori 
analysis  from  fully  resolved  DNS  at,  say,  twice  the  grid  resolution  of  the 
LES.  For  specificity,  presume  that  both  methods  exploit  Fourier  spectral 
methods.  The  comparison,  of  course,  should  occur  at  the  same  physical 
cut-off  in  wavenumber  space  (Fig.  1).  Thus  (^c)DNS  ~ (^c)LES'  however, 
because  (Ax*)kgg  = 2(Ax*)pj^g,  for  the  same  physical  cut-off,  the  dimen- 
sionless cut-offs  must  have  the  ratio  (Cc)dNS/(£c)lES  = Without  a 
tunable  filter,  the  comparison  cannot  be  made  at  the  same  physical  cut-off. 

Third,  Ad  hoc  attempts  at  tuning  the  three-point  top-hat  filter,  for 
example,  can  have  unintended  consequences.  By  naively  extending  the  filter 
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width  A*  from  2A.t*  to  4Ax*,  the  transfer  function  changes  from  low-pass 
to  U-shaped  (Fig.  5).  It  is  therefore  highly  preferable  to  begin  with  a tunable 
low-pass  filter  whose  behavior  is  suitable  at  all  values  of  dissipation  (cut- 
off). The  Pade  filter  shown  in  Fig.  6,  for  example,  is  such  a filter. 

To  summarize,  in  LES,  the  dimensional  cut-off  k*  is  determined  by  phys- 
ical considerations,  but  the  grid  resolution  Ax*  is  mandated  by  numerical 
ones.  Tunable  (one-parameter)  filters  permit  the  cut-off  and  resolution  to 
be  specified  independently,  according  to  their  respective  criteria. 

Conclusion 

Most  SGS  models  for  LES  have  been  developed  on  the  basis  of  physical 
rather  than  mathematical  considerations.  However,  it  is  shown  that  the  ex- 
act residual  stress  is  completely  determined  by  the  mathematical  properties 
of  the  grid  filter;  that  is,  by  its  order,  its  cut-off  wavenumber,  and  the  shape 
of  its  transfer  function. 
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Figure  5.  Transfer  functions  of  three-point 
top-hat  filters. 
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Figure  6.  Transfer  functions  of 

one-parameter  family  of  Pade  filters. 


Figure  4-  Planar  contours  of  exact  (E)  and 
modeled  (Ml-Smagorinsky,  M2-similarity, 
M3-gradient)  SGS  stresses  (m  compo- 
nent). 


